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ABSTRACT 



Computer simulation studies of two frequency domain adaptive beamforming algo- 
rithms are presented. The two algorithms are the frequency domain adaptive 
beamforming modified least-mean-square algorithm, and the proposed new frequency 
domain adaptive beamforming algorithm based on modified adaptive linear prediction- 
error filtering. The simulation studies were conducted to determine the multiple 
broadband target localization capability and the full angular coverage capability of the 
two algorithms. The number of iterations that the adaptive algorithms took to reach a 
minimum estimation error was determined. The algorithms were evaluated at several 
signal-to-noise ratios. Finally, using the results of the simulation studies a comparison 
between the two algorithms is performed. 
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I. INTRODUCTION 



One of the major research items in underwater acoustic signal processing is to de- 
velop new sonar signal processing algorithms capable of quickly and accurately solving 
target localization problems. Present technology dictates that several lines of bearing to 
a target must be obtained before a sonar fire control solution can be computed. Ob- 
taining these lines of bearing is a time consuming and often dangerous task due to the 
increased probability of counterdetection and, as a result, evasive maneuvering and de- 
fensive action on the part of the target. A sonar system capable of providing timely, 
accurate target localization while minimizing own ship maneuvering would result in 
longer firing ranges and, therefore, a reduction of the threat to one's own ship. 

Frequency domain beamforming is accomplished by applying appropriate phase 
shifts at the sensor outputs of an array to account for the relative propagation delays 
of a signal from a particular direction. The phase-shifted signals from all sensors are then 
added together coherently to realize the full array gain. Discrete Fourier Transform 
(DFT) beamforming is the usual method of determining the direction of arrival of a 
plane wave signal. A discrete number of direction cosine bins are formed and each bin 
corresponds to a discrete direction. If the number of direction cosine bins is large, very 
fine spatial resolution can be obtained. The phase shifts needed to cancel the relative 
propagation delays can be determined adaptively. 

The research performed in this thesis is a continuation of the work of Ziomek and 
Belirle [Ref. 1] and fully evaluates the capabilities of the complex, modified least-mean 
square, frequency domain, adaptive beamforming algorithm that they developed via 
computer simulation studies. Also, a new frequency domain adaptive beamforming al- 
gorithm based on modified adaptive linear prediction-error filtering is proposed, which 
promises better resolution in low signal-to-noise ratio environments especially when 
multiple targets share common spectral lines [Ref. 2,3]. The computer simulation studies 
were designed to test the two algorithms specifically for their multiple broadband target 
localization capability, their full angular coverage capability, and their angular resol- 
ution as a function of the input signal-to-noise ratio (SNR) at a single element in the 
array. Using lower SNR values than those used by Ziomek and Behrle [Ref. 1] in the 
evaluation of the modified LMS algorithm, and changing some assumptions, such as the 
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number of array elements, sampling rate, and number of iterations, a comparison was 
performed between the two algorithms. 

Each target was modeled as a broadband sound source. As a result, the frequency 
spectrum of the output signal from each element of the planar array contains several 
frequency components. An estimate of the bearing and depression angles associated with 
each frequency component is provided as a result of processing the output frequency 
spectrum from each element in the array through the frequency domain adaptive 
beamforming algorithms. Therefore, if each target exhibits at least one unique frequency 
component (or spectral line), then all targets can be located. 

Full angular coverage is the ability to localize a target regardless of its relative po- 
sition with respect to the array. The broadside case is the easiest since it is at this posi- 
tion that the far-field beam pattern beamwidth is its narrowest. The endfire case is the 
most difficult since the far-field beam pattern beamwidth is the broadest at this point. 
The full angular coverage and multiple broadband target capabilities were tested simul- 
taneously. 

Baseline results were the first assembled. They are defined as the bearing and de- 
pression angle estimation errors (measured in degrees) as a function of harmonic num- 
ber, and the number of iterations of the two algorithms, for the "no noise" case. 
Identical cases were then run using additive, zero mean, white, gaussian noise to corrupt 
the output signal from each element of the planar array. Average bearing and depression 
angle estimation errors were plotted as a function of the input SNR at a single element 
of the planar array, and the harmonic number. 

Chapter II describes the theory used in the development of the two frequency do- 
main adaptive beamforming algorithms. First the already known modified complex 
least-mean-square (LMS) adaptive algorithm is summarized. Then the new algorithm 
which is based on modified complex adaptive linear prediction-error filtering is discussed. 
For both algorithms the construction of the direction cosine estimates, and the angle 
estimates will be presented in detail. 

Chapter III contains computer simulation results and an explanation of these re- 
sults. The results are presented as plots of the average estimation errors of the bearing 
and depression angles at three distinct signal-to-noise ratio levels (-3dB, -6dB, -9dB) for 
five different cases for each of the two algorithms. The different test cases include targets 
evaluated at broadside, endfire, random placement, and targets which share a spectral 
line. Conclusions concerning the effect of harmonic number, sampling rate, number of 
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array elements, number of iterations used in the two algorithms, and SNR are made. The 
Appendix contains tabular numerical data for all results. 

Chapter IV will present final conclusions and recommendations for future research. 
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II. ANALYSIS 



This chapter will present a brief analysis of the two frequency domain adaptive 
beamforming algorithms for planar arrays which are evaluated in this thesis. The goal 
of both algorithms is to localize multiple broadband targets by processing the output 
electrical signals from a planar array of sensors to provide estimates of direction (both 
bearing and depression angles) and frequency content of the multiple acoustic fields in- 
cident upon the array. The development of the modified complex least-mean-square 
adaptive beamforming algorithm for planar arrays is a summary of the analysis section 
of a paper by Ziomek and Behrle [Ref. 1]. However the modified complex adaptive lin- 
ear prediction-error filter beamforming algorithm for planar arrays is presented for the 
first time. As a theoretical background for both algorithms, the principles of adaptive 
beamforming for planar arrays in the frequency domain are presented in the next sec- 
tion. 



A. FREQUENCY DOMAIN ADAPTIVE BEAMFORMING FOR PLANAR 
ARRAYS 

Consider a planar array of M x X (odd), equally spaced, point source elements lying 
in the XY plane where M and N are the total odd number of elements in the X and Y 
directions, respectively. Let the random output electrical signal received at time instant 
/ and element (m,n) in the array be given by 



r{Ltn,n) = y(l,m,n ) + l = -Z/,...,0,...,Z/ 

m = — A/' ,...,0,...A/' (2.1) 

n = — A r ',... f 0,...iY' , 

where y(l,m,n) and n(l,m,n ) are the deterministic signal and random receiver noise, 
respectively, 

L' = {L-\)12 (2.2) 
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M' = (M - l)/2 



(2.3) 



i\’ = {N- l)/2 



(2.4) 



and 



L>2K+l (2.5) 

is the total number of time samples that must be taken per element in order to avoid 
aliasing when the signal y(l,m,n) is composed of K harmonics [Ref. 4]. If we designate 
the length of the data record (i.e., the fundamental period ) recorded at each element in 
the array as T 0 seconds, then the fundamental frequency (i.e., the FFT bin spacing ) is 
f 0 = -y- Hz, the highest frequency component that can be calculated in y(l,m,n) is Ky^ 

f 0 

Hz, and the output received signals r(l,m,n ) must be sampled at a rate of 



f 5 = L/T 0 samples, sec, (2.6) 

where L is given by Eq. (2.5) [Ref. 4). Taking the DFT of Eq. (2.1) with respect to the 
time index / yields the complex frequency domain samples. 



R(q,m,n ) = Y{q.m,n ) + N(q,m,n), q = — L',...,0,...,Z/ 

m = (2.7) 

n = -,V',...,0,...iV , 

where the index q represents the harmonic number. If the acoustic field incident upon 

V * YX 

the planar arrray is a single general plane wave g (t + - c — ) propagating in the ± n 0 di- 
rection (see Fig. 1), where « 0 is a unit vector and g(t) is an arbitrary baseband function, 
then it can be shown that the output electrical signal at time instant / and element (m,n) 
in the array is given by [Ref. 4 p. 160]: 



„ , u 0 md x + v () nd, 

y{Lm,n) = g{!T s + ) . 



( 2 . 8 ) 
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, ) 

Figure 1. General plane-nave field g (/ + — ) propagating in the + direction 
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The corresponding frequency spectrum is given by [Ref. 2] 



,_j2nqf 0 u 0 md x j2nqf 0 v 0 nd y 

} (q,m,n) = Lc q exp( + ~ c ) exp( + ) 



(2.9) 



where 



c q = Y(q,0,0)IL, q = — L',...,0,...,L' (2.10) 

are the complex Fourier series coefficients that can be used to represent g(t) by a finite 
Fourier series with K harmonics during the time interval |r| < — ~,f 0 = -4r- , 

l 1 q 



w 0 = sin B 0 cos \j/ 0 



( 2 . 11 ) 



and 



v 0 = sin B 0 sin \p 0 (2. 1 2) 

are the dimensionless direction cosines with respect to the X and V axes, respectively; 
B 0 and rp „ are the depression and bearing angles, respectively; d„ and d y are the interele- 
ment spacings in meters in the X and Y directions, respectively; and c is the speed of 
sound in meters per second. 

Equation (2.9) is the basis for the frequency domain adaptive beamforming algo- 
rithms to be presented in this thesis. The algorithms presented use a planar array to 
estimate both B 0 and for each harmonic q present in the multiple incindent plane-wave 
fields via frequency domain adaptive beamforming. Therefore, we shall use the following 
generalized version of Eq. (2.S) 



y(l,m,n) = ^g k (lT s 



_ u 0k md x + v ok nd y 

j + c ) 



(2.13) 



which represents the output electrical signal at time instant / and element (m,n) in the 
planar array due to several incident plane-wave fields, where g k (t), w 0A , and % are the ar- 
bitrary time function and direction cosines associated with the kth sound source (target), 
respectively. 
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B. FREQUENCY DOMAIN MODIFIED COMPLEX LMS ADAPTIVE 
BEAMFORMING ALGORITHM 



The frequency domain adaptive beamforming algorithm to be discussed is based 
on processing the output complex frequency domain data R(q,m,n) from all M xN ele- 
ments in the planar array. We begin by defining the complex estimation error as 



e(q) = s(q) - s(q) 



(2.14) 



where 



\r y 

*(<?)= J ^ y] \R{q,m.n)\ e\pl+jLR{q.0.0)'] 

m=—M' n=—.V 



(2.15) 



is the reference signal; 



a v y 

s(q) = T \/V 2 ^ c{q.m)d(q,n)R{q,m,n) 

m=-\r n=-Y 

c T (q)R(q)d(q) 

LMX 



(2.16) 



is the estimate of s(q); c(q.m) and d(q,n) are the unit magnitude complex weights in the 
X and Y directions, respectively; £(</) and d(q) are the M x 1 and N x 1 complex weight 
vectors in the X and Y directions, respectively; given by 



dq) = lc(q,-M'), 



,c(q.0),...,c(qMY 



(2.17) 



and 



d{q) = [ d(q , - X%..,d(^),...,d(q,y)l T 



(2.18) 
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respectively; and 



M<?) = 



R{q, - A/', - X')...R(q, - M',0)...R{q, - A/', A") 
R(qA-X’) - R(t 7.0,0) ... £(</,<), A") 

R(q,\i’, — A"') ... R{q,M‘ ,0) ... R(q,M',iY) 



(2.19) 



is the MxN complex data matrix. Next, define the (M + X) x 1 complex weight vector 
vv(^) as follows: 






dq) 

diq) 



( 2 . 20 ) 



Therefore. 



c{q) = A w{q) 



( 2 . 21 ) 



where 



A=l 1 IQ] 

A/.vA/ A/x.V 



( 2 . 22 ) 



is a M x (M + X) matrix, / is a M x M identity matrix, and 0 is a M x X null matrix 
as indicated; and 



d{q) = B w{q) 



(2.23) 



where 



Q I I 

NxM NxN 



(2.24) 



is a X x (M + X) matrix, 0 is a X x M null matrix, and / is a X x X identity matrix as 
indicated. Substituting Eqs. (2.21) and (2.23) into Eq. (2.16) yields 
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(2.25) 



a f , w T {q)Z(q)w{q) 

S[q) = UIX 

where 

Z{q) = A T R{q)B (2.26) 

is a (M + N) x (M + N) complex matrix. The complex weight vector that minimizes the 
mean-square error £{|e(g)| 2 } is given by [Ref. 1] 

w /+1 (^) = vv / (?) + 2 ¥ / ? )[Z( ? ) + Z r ( ? )]V, / =0,1,2,.... (2.27) 

where 

e i(q) = s(q) - s£q) ■ (2.28) 

is the estimation error after the ith iteration, 



m = 



Ki(q)Z(qWq) 

LMX 



(2.29) 



is the estimate of s(q) after the ith iteration, and, the step-size parameter /t, is given by 

M/ = Mo = + ^r 1 . ‘ = 0,1,2, (2.30) 

where ju 0 is constant and a] and a 2 n are the signal and noise power, respectively, at the 
center element in the array. After each iteration, each component of the complex weight 
vector w,-j(< 7) is normalized by its respective magnitude in order to maintain unit mag- 
nitude. 

Once the complex weight vector vv,.,(^) converges to a steady-state value Ws S (q), the 
steady-state complex weight vectors £s S (q) a ^d d ss {q ) can be obtained from Eqs. (2.21) 
and (2.23), respectively. The estimates of 6 0 (q) and i p 0 (q) at each harmonic q are given 
by [Ref. 3] 
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(2.31) 



Uq ) = sin -'[{[^)] 2 + [ tf 5 (?)] 2 }' /2 ], q * 0 



and 



<Ao(?) = tan 



-1 



a LS 



(?) 



a LS / \ 

«o (?) 



q ¥* 0 



(2.32) 



where ut s (q) and ii s (q) represent estimates of the direction cosines obtained by using a 
least-squares fit to the "unwrapped" steady-state phase weights d l ss (q,m) and 4>% s (q,n), 
respectively. Note that in the absence of noise 



Qss(q- m ) = ± 



M(q)nui x 

c 



q = -L’,...,0,...,L’ 



(2.33) 



and 



</>ss(?-«) = ± 



o Y o(q) nd x 

c 



q = — L',...,0,...,L' 
n = — .V,....0 A 7 



(2.34) 



The steady-state phase weights need to be "unwrapped", that is, allowed to take on val- 
ues outside the closed interval [ — n , «], in order to ensure full angular coverage (i.e., 
0 < 6 0 (q) < — and 0 < \Jt 0 {q) < 2n) and correct depression and bearing angle estimates, 

A A ~ 

6 0 (q) and ip 0 (q), respectively. Finally, substituting Eqs. (2.15) and (2.16) into Eq. (2.14) 
yields the following expression for the steady-state estimation error: 



\r ,v 

e$s(?)- LV/V Z £ \R(q,m,n)\ exp[ +jLR(q, 0,0)] 
m=-\r n=-y 

- exp{ +jlds S {q,m) + 4>ss(q,n) + LR{q,m,n)']} 



(2.35) 



where 
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R{q,m,n) = \R{q,m,n)\ exp[ +jLR(q,m,n)'] 



(2.36) 



II” W 

= «55(<7^) ex p[ +/0 sj(< 7^)] = ex PC +jdss(q> m )'] ( 2 .37) 



and 



«ss(<7> w ) = 6ss(?> w ) ex PC +^55(<?,«)] = ex P[ +j4>ss(q ,")] 



(2.38) 



where a ss (q,m)= 1 and b ss (q,n ) = 1 are real, unit magnitude, amplitude weights and 
05 s(< 7 ,m) and 4>f s (q,n) are real, "wrapped" phase weights. The procedure for the unwrap- 
ping of the steady-state phase weights in order to obtain least-squares estimates of the 
direction cosines at each harmonic is described by Behrle [Ref. 5]. 

C. FREQUENCY DOMAIN MODIFIED COMPLEX ADAPTIVE LINEAR 
PREDICTION-ERRROR FILTER BEAMFORMING ALGORITHM 

Consider a planar array of M x N (odd), equally spaced, point source elements lying 
in the XV plane where M and N are the total odd number of elements in the X and Y 
directions, respectively. Figure 2 shows a linear prediction-error filtering structure in the 
frequency domain (i.e., after the operation of a time domain DFT). The frequency do- 
main adaptive algorithm to be presented is based on processing the output complex 
frequency domain data R(q,m,n) given by Eq. (2-7). Since the target localization prob- 
lem requires only phase informations in our algorithm, we processes phase only data. 

Let the one-step forward linear prediction errors for harmonic number q in the X 
and Y directions, be given by 



= s x (q,i,n) - s x {q,i,n), n = - A\...,0,...,A" 



(2.39) 



and 
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1 i i + l i+2 i+3 i+P M x 




Figure 2. Frequency domain linear prediction-error filter structure 

respectively, with 

A T ' = ~ V ~ 1 (2.41) 

and 

Af = — f — , (2.42) 

where i is the number of the reference element in the direction of processing; • ) and 
5 ,( •) are the reference signal and the estimate of the reference signal in the X direction, 
respectively; and • ) and £,( • ) are the reference signal and the estimate of the reference 
signal in the Y direction, respectively. Each one of the above signals is defined as follows: 



s x (qj,n) = R(q,i,n), 



(2.43) 



p 

a 1 V 1 J t 

s x (q,i,n) = ~p 2^h x (q,p,n)R(q,i + p,n) = — h x (q,n)L(q,i,n), 

p= i 



(2.44) 
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s y (q,m,i) = R(q,m,i), 



(2.45) 



and 



p 

A 1 \ ^ J j' 

s y (q,m,i) = — 2j^W)R{qW + p) = ~pt y {q,m)r{q,m,i) (2.46) 

P = i 

where h x {q.n) and h y (q,m) are the P x 1 complex weight vectors in the X and Y directions, 
respectively , and r{q,i,n ) and r(q,m,i ) are the P x 1 complex frequency domain data vec- 
tors in the X and Y directions, respectively, given by 



h x {q.n) = [h x {q,l,n), h x (q,2,n), ,h x {q,P,ri)Y, (2.47) 

liq.i.n) = [R{q,i + \,n),R{q.i + 2 .«), ,R{q , / + P,n)Y , (2.4S) 

ky{q.m ) - C hy(q.inA ). h y (q,m,2) ,h y (q,m,Py] T , (2.49) 

and 

r{q,m.i) = \_R(q,m,i + 1 ),R(q,m,i + 2) ,R(q, m,i + P)] r (2.50) 

where P is the order of the linear prediction-error filter. 

Substituting Eqs. (2.43) and (2.44) into Eq. (2-39) and, substituting Eqs. (2.45) and 
(2.46) into Eq. (2.40) yields 



p 

F Px {q,i,n) = R(q,i,n ) - y ^ h x {q,p,n)R{q,i + p,n), n = - A%..,0,...,A r ', 

P=1 



(2.51) 



and 
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m — — V/'. (2.52) 



p 

F Py (q, m ,i) = 0 - -jr ^‘ y (q,m,p)R(q,m,i + p), 

p= 1 



Since we are working with phase only data, Eqs. (2.51) and (2.52) become: 



= exp[ +jLR{q,i,ny] 

P 

- -J ^Sjh x {q,p,n)\ exp[ +jLh x (q,p,n)'] exp[ +jLR{q,i + />,«)], (2.53) 

p= i 

n = — A r ',...,0,...,iV' 



and 



Fp(qqn.i) = exp[ t/Z_ R{q, m,/)] 

P 

- y ^Jfy(?,m./7)| exp[ +jLh y {q,m,p )] exp[ +jLR{q,m,i + /?)]. (2.54) 

P =\ 

m = -3f ,...,0,...,3/'. 

From Eqs. (2.53) and (2.54) it can be seen that if each component of the complex 
weight vectors have unit magnitude, then the forward linear prediction error become 
zero when 



exp[ +jLR{qJ,n)2 = y V exp[ +jLh x {q,p,n)~] exp[ +jLR{q,i + /?.«)], 



p=\ 



(2.55) 



n = — A',...,0,...,iV' 



and 



r 

exp[ +jLR(q.m,i)2 = -4- ^ exp[ +jLh y (q,m,p)] exp[ +jLR{q,m,i + />)], 



p = l 



(2.56) 



am - — 3f ,...,0 ,...,3f . 
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Equations (2.55) and (2.56) will be valid only if the phase weights cophase the data at 
the output of each element in the prediction-error filters. 

The algorithm used for the optimum estimation of the complex weight vectors is a 
modified, complex least mean-square (LMS) algorithm. The estimates of the complex 
weight vectors in the X and V directions after the k-th iteration in X and Y directions 
is given by 



\_h x {q.n)'] k+l = ih x (q.n)'] k + 2/x A [F Pj[ (?, /»]*[£*(?, /.«)], 

n — — iV',...,0,...A’' (2-57) 

A — 0, 1,2,3,.... 

and 



= lh y {q.m)'] k + 2 ^[/>^.m,/)] A [r*(^.m,/)], 

m = - d /',....0 M' ( 2 - 58 ) 

A- = 0.1, ir- 



respectively. where the index k denotes the k-th step of iteration, F Px ( • )* and. F t { • )* are 
the forward linear prediction errors in the X and Y directions, respectively, after the k-th 
iteration: £*( • ) is the complex conjugate of the data vector; and the step-size parameter 
,u k is given by 



Mfc = Mo = (°y + °n) \ A = 0,1,2,... (2.59) 

where is a constant, and it is equal to the inverse of the sum of the signal and noise 
power, a] and a 2 n , respectively, at the center element of the array. After each iteration 
of the algorithm, each component of the complex weight vectors is normalized by its 
respective magnitude in order to maintain unit magnitude. 

Since we are working with multiple linear arrays to obtain estimates of the targets' 
locations, we use average, steady-state, complex weight vectors in the X and Y di- 
rections, that is, 
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(2.60) 



v 

h x Jq) = jf V/^,/ 7 ) 

/7=-.V 

and 




(2.61) 



Note that each of the P x 1 average, steady-state, complex weight vectors in the X and 
Y directions is given by 



and 



= l)AJ<7-2), ,h.{q,P)y 



(2.62) 



4’J<7) = I XJ?- 1 )’ VJ?- 2 )’- 



AJq.P)?. 



(2.63) 



Next, if we now use the average, steady-state, complex weight vectors in the X and 
Y directions, then the one-step, steady-state, forward linear prediction errors for har- 
monic number q in the X and Y directions, are given by 



r 

(<?,'» = R{q,i,n ) - — ^h Xovg (q,p)R(q,i + p,n), n = - N’ , (2.64) 

p= i 



and 

p 

F Pg {q,m,i ) = R{q,m,i ) - -y '^Ji yfn i,q,m,p)R{q,m,i + p), m = - M’ ,...,0,...,Af . (2.65) 

p=i 



If we define the forward spatial Z-transform of g(i) as 
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(2.66) 



DO 

G(Z) = Z{g(i)} = £ 

i — — oo 



then 



Z{g{i + p)) = G{z)z ~ p , (2.67) 

and, upon taking the forward spatial Z-transform of Eqs. (2.64) and (2.65) with respect 
to i, we obtain 



r 

F Px (q.z x .n) = R(q.z x ,n)ll - j? Y h ^'P) z 7^ « = - AV-O,...,*', (2.68) 



P = 1 



and 



f 

F Pji (q,m,z y ) = R(q,m,Zy)l 1 - y m = “ M' ■ (2.69) 

p= 



From Eqs. (2.6S) and (2.69), the transfer functions of the one-st6p forward linear 
prediction-error filters (analysis filters) in the X and V directions are given by 



HpE y ( ( i- : x’ n ) — ' 



F p (q.z x ji) 



R{q,z x .n) 



p 

”~~ iV 0 (2 - 70) 

p = i 



and 



F p (q, m ,Zy) | 

UpE^RMy) = = 1 - “p YhjA'ti 2 * ' m = ~ ( 2 - 71 ) 



Since the average complex weights in the X and Y directions are independent of the in- 
dices n and m, the transfer functions are also independent of n and m, that is, 
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p 

H PEx {q,z x ) =l-y Vh.Jq.p):^, (2.72) 

P= 1 

and 

p 

Hpe,^ - 1 " T 2X,(«’) 2 7' < 2 ' 73 ) 

p=l 



Therefore the transfer functions of the inverse filters (synthesis filters) in the X and Y 
directions are given by 



and 



= r 1 (2-74) 

1 - 77 

P= 1 

««,(«•*,) . (2.75) 

P= 1 



In order to obtain the spatial-frequency response (angular response) of the filter 
described by Eq. (2.74), let 



+j2nud x 

z x = exp( +j2nf x d x ) = exp( — - ). 



(2.76) 



Substituting Eq.(2.76) into Eq.(2.74) yields 



HpeXQ’ 14 ) 



1 - y exp( 



—j2npud x 



p = 1 



(2.77) 
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Similarly, if we let 



Zy = exp( +j2nfyd y ) = exp( 



+j2nvd y 



(2.78) 



in Eq. (2.75), we obtain 



Hply 7.V) = 




(2.79) 



Equations (2-77) and (2-79) are then evaluated as functions of the direction cosines 
u and v, respectively, for a given harmonic q. The values of u and v that maximize Eqs. 
(2-77) and (2-79), respectively, then become the linear prediction-error estimates u£ E {q) 
and \i E {q) of u 0 (q) and v 0 (c/) . The estimates of 6 a (q) and <j/ a {q) at each harmonic q are 
given by [Ref. 3] 



where ii^ E (q) and v^ E (q) are the estimates of the direction cosines obtained by using the 
adaptive linear prediction-error algorithm. 




(2.80) 



and 




q ¥= 0 



(2.81) 
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III. SIMULATION RESULTS 



Computer simulation results for five test cases are presented in this chapter. The test 
cases were designed to test both of the algorithms multiple broadband target localization 
capability, full angular coverage capabilities, and their angular resolution as a function 
of the input SNR at a single element in the array, harmonic number, and the number 
of iterations. 

For simplification purposes, in all of the simulation results we use the following 
notation: 

• L.MS: When the frequency domain modified complex Least-Mean-Square adaptive 
beamforming algorithm is used. 



• ALP: When the frequency domain modified complex adaptive linear prediction- 
error filtering beamforming algorithm is used. 

The test cases are described below: 

1. Case 1 - A single broadband target radiating six harmonics, located at broadside 
relative to the planar array; 

2. Case 2 - Three broadband targets, each radiating two unique harmonics, located 
at random positions; 

3. Case 3 - A single broadband target radiating six harmonics, located at endfire rel- 
ative to the planar array; 

4. Case 4 - Three broadband targets. Two targets are located in the same plane and 
have a common harmonic. Three subcases were examined. Each subcase used a 
different harmonic as the common spectral line. 

5. Case 5 - Three broadband targets, located at random positions. Two targets have 
a common harmonic. Three subcases were examined. Each subcase used a different 
harmonic as the common spectral line. 

The simulation results of these five cases are based upon processing the output 
electrical signals from a 1 1 x 1 1 planar array of equally spaced hydrophones. The 
acoustic field incident upon the planar array was, in general, the sum of several plane- 
wave fields travelling in different directions. Each plane-wave field consisted of an arbi- 
trary number of harmonics (spectral lines) emanating from one of the broadband sound 
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sources. As a result, the output electrical signal at each element of the planar array was 
composed of an arbitrary number K of harmonics, all with identical amplitudes of unity. 

The fundamental frequency for all test cases was chosen to be 7 = 1000 Hz. Via Eq. 
(2-6), the fundamental period at each element of the array was 1 millisecond. The sam- 
pling parameter, S, was set equal to 10. The number of time samples taken per element 
of the array is given by 



L = SK+ 5 (3.1) 

where L is the total number of samples taken per element, K is the total number of 
harmonics present in the signal, and S is the sampling parameter. With S= 10 and K = 6, 
L = 65 time samples per element. 

In the computer simulations presented in this thesis, the interelement spacing is 



d x dy 



(3.2) 



where the minimum wavelength is 



'•min = ~T ~ . (3.3) 

J max 

and 

/ max =A7 0 Hz (3.4) 

Substituting Eqs. (3.3) and (3.4) into Eq. (3.2) yields 




(3.5) 



For the ALP algorithm, the order of the forward prediction-error filter was 10. The 
reference element was the first element of each linear array and the number of linear 
arrays that have been processed in both directions was 11. Initial simulation tests indi- 
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cated that the number of iterations required for the LMS algorithm to approach a 
steady-state error was 100, while for the ALP algorithm 1000 iterations were required. 

Baseline, or "no noise" test case results were the first generated to ensure that both 
algorithms were working properly in a noise-free environment. Figure 3 depicts the 
noise-free, time domain, received signal at element (1.1) for Case 2. For each baseline 
test case, bearing and depression angle estimation errors (measured in degrees) were 
obtained by running the computer simulation once allowing 100 iterations for the LMS 
algorithm and 1000 iterations for the ALP algorithm. 

Following compilation of the baseline results for the first three cases, identical test 
cases were run using additive, wide-sense stationary, zero mean, white, gaussian noise 
samples to corrupt the time samples of the received signal. Figure 4 depicts the time 
domain, received signal at element (1.1) for case 2 for a SXR= -9dB. For each test case, 
and for a given input SNR at a single element of the array, average bearing and de- 
pression angle estimation errors were obtained by running the computer simulation 10 
times. During each run. the LMS algorithm was allowed 100 iterations whereas the ALP 
algorithm was allowed 1000 iterations. 
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V TIME INSTANT L 

CASE: LMS2 

INPUT SNR AT ELEMENT (1,1): NO NOISE 



Figure 3. Real Received Signal at Element (1,1), for Case LMS2 : No Noise 
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CASE: LMS2 

INPUT SNR AT ELEMENT (1,1): -9.0 DB 



Figure 4. Real Received Signal at Element (1,1), for Case LMS2 : SNR = -9dB 
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A. CASE 1 



Case 1 placed a single broadband target at broadside relative to the planar array 
(i.e., 6 0 = 0°). This was considered the simplest case for both of the algorithms since it 
is at broadside that the far-field beam pattern has its narrowest beamwidth and, as a 
result, both algorithms should provide good angular resolution [Ref. 4]. The general 
plane-wave field radiated by the target consisted of six harmonics. 

Figures 5, 6 and 7 present the average estimation errors of the depression angle 
( e 0g ) versus SNR for Case LMS1 at 1000 Hz, 3000 Hz and 6000 Hz, respectively. Results 
for bearing angle estimation errors in the broadside case are irrelevant since the target 
is directly above the array, and the bearing angle has no meaning. From these figures, 
we can see two trends. Increasing the SNR results in smaller magnitudes of the average 
estimation errors of the depression angle and as the harmonic number(q) increases, the 
magnitude of the estimation error decreases. These two trends were expected. As we 
increase the SNR, the noise components of the received signals are less dominant. An 
increase of the harmonic number (q) represents an increase in the frequency and, as a 
result, the beamwidth of the far-field beam pattern decreases. The decreased beamwidth 
increases the angular resolution of the algorithm and, as result, decreases the estimation 
error. 

Figures 8, 9 and 10 present estimates of the direction cosines U 0 and V 0 for Case 
ALP1 at SNR=-3 dB and frequencies 1000 Hz, 3000 Hz, and 6000 Hz, respectively. 
Similarly. Figures 11,12 and 13, and 14, 15, and 16 present estimates of the direction 
cosines C 0 and V 0 for Case ALP1 with the same order of frequency as above and for 
SNR=-6 dB and SNR = -9 dB. respectively. As we discussed in Chapter 2, the location 
of the peaks of the magnitude of the synthesis filter transfer function in dB give the U 0 
and V 0 estimates. From these figures we can clearly observe that as the SNR decreases, 
the difference between the peak and the background noise in the transfer function plots 
decreases. Also, the harmonic number is important, because increasing the harmonic 
number produces sharper peaks in the transfer function plots resulting in more accurate 
estimation. 

Figures 17, 18 and 19 present the average estimation error of the depression angle 
( e 6o ) versus SNR for Case ALP1 at 1000 Hz, 3000 Hz and 6000 Hz. The trends which 
were present in Case LMS1 are apparent here, that is, as SNR increases, the magnitude 
of the estimation error decreases and, as the harmonic number (q) increases, the mag- 
nitude of the estimation error decreases. 
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Finally Tables 1 and 2 in the Appendix A present the numerical data for Cases 
LMS1 and ALP1 respectively for all six harmonics. Comparing the two algorithms us- 
ing these tables, we can conclude that for Case 1, both algorithms give small magnitude 
estimation errors in low SNR with a little better performance for the LMS1 particularly 
at 1000 Hz. 
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CASE: LMS1 

FREQ: 1000.0 HZ M- 1 1 N- 11 NITER- 100 



' t 1 

Figure 5. Average Estimation Error of Depression Angle vs. SNR at 1000 Hz for 
Case LMS1. 
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AVERAGE ESTIMATION ERROR (DEG) 
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Figure 
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- 9.0 - 6.0 - 3.0 

INPUT SNR (DO) RT R SINGLE RRRRY ELEMENT 

CRSEtLMSl 

FREQ: 3000.0 HZ M- 1 1 N- 1 1 NITER- 100 



6. Average Estimation Error of Depression Angle vs. SNR at 3000 Hz for 
Case LMS1. 
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AVERAGE ESTIMATION ERROR (DEG) 

- 10.0 - 8.0 - 6.0 - 4.0 - 2.0 0.0 2.0 4.0 6.0 8.0 10.0 




CASE: LMS1 

FREQ: 6000.0 HZ M- II N- 1 1 NITER- 100 



Figure 7. Average Estimation Error of Depression Angle vs. SNR at 6000 Hz for 
Case LMS1. 
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Figure 8. Estimates of Direction Cosines UO and VO, Case ALP1: SNR -3dB and 
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CASE J flLPl SNR- -3.0 00 

FRE0» 3000.0 HZ h- 11 N- II I- 1 P-10 NITER- 1000 



Figure 9. Estimates of Direction Cosines U0 and VO, Case ALP1: SNR = -3dB and 
q = 3 
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Figure 10. Estimates of Direction Cosines UO and VO, Case ALP1: SNR = -3dB and 

q = 6 
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Figure 11. Estimates of Direction Cosines UO anti VO, Case ALP1: SNR = -6dB and 
q= 1 



34 



MAGNITUDE IN 08 MAGNITUDE IN DB 





















I\ 






































\ 




































I 


1 




















y 












\ 




/ 

A 




\ 


s, 






\ 






















\ 




/ 












































































































































































































- 









- 1 . 0 - 0.9 - 6.8 ' 0.7 - 0.6 - 0 .S - 0.4 - 0.3 - 0.2 ' 0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 o ‘.7 0.8 0.9 1.0 

U 



















i 


\ 


































1 


\ 


































in 


\ 


























/ 








/ 


\ 


























/ 




\ 




/ 


\ 




/ 






— 





























































































































































































' - 1 . 0 - 0.9 ' 0.8 - 0.7 - 0.8 - 0.5 - 0.4 - 0.3 - 0.2 ' 0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



CASE: RIP 1 SNR- -6.0 DB 

FREOi 3000.0 HZ H- 1 1 N- II 1-1 P- 10 NITER- lOOu 

Figure 12. Estimates of Direction Cosines UO and VO, Case ALP1: SNR = -6dB and 
q=3 
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Figure 13. Estimates of Direction Cosines UO and VO, Case ALP1: SNR = -6dB and 
q=6 
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Figure 14. Estimates of Direction Cosines UO and VO, Case ALP1: SNR = -9dB and 
q= I 
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Figure 15. Estimates of Direction Cosines UO and VO, Case ALP1: SNR = -9dB and 
q=3 
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CASE: ALP 1 SNR- -9.0 00 
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Figure 16. Estimates of Direction Cosines UO and VO, Case ALP1: SNR = -9dB and 



q = 6 
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Figure 17. Average Estimation Error of Depression Angle vs. SNR at 1000 Hz for 
Case ALP1. 
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Figure 18. Average Estimation Error of Depression Angle vs. SNR at 3000 Hz for 
Case ALP1. 
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B. CASE 2 



Case 2 placed three broadband targets at random locations. The general plane-wave 
field radiated by each of the targets contained two unique harmonics . Target 1 was lo- 
cated at ( 6 0 = 49°, \{/ 0 = 38°) and radiated harmonics 1 and 6. Target 2 was located at 
(0 O = 5°, \j/ 0 = 137°) and radiated harmonics 2 and 5. Target 3 was located at (9 0 = 11°, 
\p o = 307°) and radiated harmonics 3 and 4. Since there are three incident plane-wave 
fields, each containing two unique harmonics, the output electrical signal from each el- 
ement in the array exhibits a total of six harmonics. 

Figures 20 through 25 present the average estimation errors of the depression (c<, 0 ) 
and bearing (e+ ) angles versus SNR for Case LMS2 for all the harmonics in increasing 
order of the harmonic number (q), respectively. The same trend about the SNR as in 
Case 1 holds here, that is increasing the SNR results in smaller magnitude of the average 
estimation error. These general results can be explained using the same arguments pre- 
sented for Case 1. 

Figures 26 through 43 present estimates of the direction cosines U 0 and V 0 for all the 
harmonics placed in increasing order of the harmonic number (q) at SNR's -3 dB, -6 dB 
and -9 dB, respectively. Again, the two trends observed in Case ALP1 are apparent here, 
that is, decreasing the SNR decreases the difference between the peak and the back- 
ground noise and increasing the harmonic number (q) produces sharper peaks in the 
transfer function plots. 

Figures 44 through 49 present the average estimation errors of the depression ( e So ) 
and bearing (^ ) angles versus SNR for Case ALP2 for all the harmonics in increasing 
order of the harmonic number (q), respectively. The same trends as in L\IS2 hold, with 
the exception that now the magnitude of the average estimation error for both angles is 
smaller than in LMS2 . Note the difference in scaling between the error plots for Case 
LMS2 and Case ALP2. 

Finally, Tables 3 and 4 in the Appendix A present the numerical data for Cases 
LMS2 and ALP2, respectively, for all six harmonics. Comparing the two algorithms 
using these tables, we can conclude that for Case 2, the ALP algorithm gives smaller 
magnitude estimation errors in most instances than the LMS algorithm especially for 
low SNR. 
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Figure 20. Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 1000 Hz for Case LMS2. 
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Figure 21. Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 2000 Hz for Case LMS2. 
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AVERAGE ESTIMATION ERROR (OEG) 
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Figure 22. Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 3000 Hz for Case LMS2. 
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COSE: LMS2 

FREO: 4000.0 HZ M- 1 1 N- 1 1 NITER- 100 



Figure 23. Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 4000 Hz for Case LMS2. 
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AVERAGE EST I MAT I ON ERROR(DEG) 
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Figure 24. Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 5000 Hz for Case LMS2. 
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Figure 25. Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 6000 Hz for Case LMS2. 
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CR5E:RLP2 SNR- -3.0 00 
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Figure 26. Estimates of Direction Cosines UO and VO, Case ALP2: SNR = -3dB and 
q= 1 
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FREOt 2000.0 HZ h- 1 1 N- 11 I- l P- 10 NITER- 1000 

Figure 27. Estimates of Direction Cosines UO and VO, Case ALP2: SNR = -3dB and 
q = 2 
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Figure 28. Estimates of Direction Cosines UO and VO, Case ALP2: SNR = -3dB and 
q=3 
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CflSE:flLP2 SNR- -3.0 00 

FREOs -1000.0 HZ M- II N- 1 1 l-l P- 10 NITER- 1000 



Figure 29. Estimates of Direction Cosines UO and VO, Case ALP2: SNR = -3dB and 
q = 4 
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Figure 30. Estimates of Direction Cosines UO and VO, Case ALP2: SNR = -3dB and 
q=5 
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Figure 31. Estimates of Direction Cosines UO and VO, Case ALP2: SNR = -3dB and 
q = 6 
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Figure 32. Estimates of Direction Cosines UO and VO, Case ALP2: SNR = -6dB and 
q= 1 
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Figure 33. Estimates of Direction Cosines UO and VO, Case ALP2: SNR = -6dB and 
q = 2 
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Figure 34. Estimates of Direction Cosines UO and VO, Case ALP2: SNR = -6dB and 
q = 3 
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Figure 3?. Estimates of Direction Cosines UO and VO, Case ALP2: SNR 
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Figure 36. Estimates of Direction Cosines UU and VO, Case ALP2: SNR = -6dB and 
q = 5 
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CASE*. ALP 2 SNR- -6.0 00 

FREQt 6000.0 HZ h- 1 1 N- 1 1 I- I P- 10 NITER- 1000 

Figure 37. Estimates of Direction Cosines UO and VO, Case ALP2: SNR = -6dB and 
q=6 
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CflSE:flLP2 SNR- -9.0 00 
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Figure 38. Estimates of Direction Cosines UO and VO, Case ALP2: SNR — -9dB and 
q=l 
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Figure 39. Estimates of Direction Cosines U0 and VO, Case ALP2: SNR = -9dB and 
q = 2 
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Figure 40. Estimates of Direction Cosines UO and VO, Case ALP2: SNR = -9dB and 
q=3 
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CASE: RLP2 SNR- -9.0 00 

FREO: <000.0 HZ M- 11 N- 1 1 I- 1 P- 10 NITER- 1000 



Figure 41. Estimates of Direction Cosines UO and VO, Case ALP2: SNR = -9dB and 
q = 4 



65 





CASE: ALP2 SNR- -9.0 D0 

FREO: 5000.0 HZ h- 1 1 N- 11 l-l P- 10 NITER- 1000 



Figure 42. Estimates of Direction Cosines UO and VO, Case ALP2: SNR = -9dB and 
q = 5 
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Figure 43. Estimates of Direction Cosines UO and VO, Case ALP2: SNR = -9dB and 
q = 6 
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Figure 44. Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 1000 Hz for Case ALP2. 
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Figure 45. Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 2000 Hz for Case ALP2. 
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Figure 46. Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 3000 Hz for Case ALP2. 
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AVERAGE ESTIMATION ERROR(DEG) 
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Figure 47. Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 4000 Hz for Case ALP2. 
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Figure 48. Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 5000 Hz for Case ALP2. 
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Figure 49. Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 6000 Hz for Case ALP2. 
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C. CASE 3 



Case 3 placed a single broadband target at endfire relative to the planar array (i.e., 
6 0 = 90°, i/^ = 90°). This case was considered the most difficult for both algorithms since 
it is at endfire that the far-field beam pattern beamwidth is broadest [Ref. 4]. The general 
plane-wave field radiated by the target consisted of six harmonics. 

Figures 50 and 51 present the average estimation errors of the depression ( e $0 ) and 
bearing (e H ) angles versus SNR for Case LMS3 at 1000 Hz and 3000 Hz, respectively. 
The estimation errors for the highest harmonic (6000 Hz) are not shown for Case LMS3 
since the magnitude of these errors were between 45 and 90 degrees. This poor per- 
formance for the highest harmonic at endfire can be explained theoretically [Ref. 6]. 

Figures 52, 53. and 54 present estimates of the direction cosines U 0 and V 0 for Case 
ALP3 at S\R=-3 dB and frequencies 1000 Hz, 3000 Hz, and 6000 Hz, respectively. 
Similarly. Figures 55, 56, and 57, and 5S, 59, and 60 present estimates of the direction 
cosines U 0 and V 0 for Case ALP3 with the same order of frequency as above and 
S\R=-6 dB and SNR=-9 dB. respectively. Here we can observe an ambiguity for the 
f' 0 estimates at 6000 Hz which can be explained by examining the exponential term of 
Eq.(2-79). This term for the highest frequency component creates integer multiplies of 
2k producing two peaks at ± 1 for the V 0 estimates at 6000 Hz. 

Figures 61. 62 and 63 present the average estimation errors of the depression ( e eQ ) 
and bearing (e^) angles versus SNR for Case ALP3 at 1000 Hz, 3000 Hz, and 6000 Hz, 
respectively. We can observe that in Case ALP3 at 6000 Hz we have acceptable esti- 
mates for depression angle (6 0 ). In Figure 63 the t/' 0 error plot is not shown since the 
magnitude of this error was between 45 and 90 degrees. 

Finally, Tables 5 and 6 in Appendix A present the numerical data for Cases LMS3 
and ALP3 respectively for all six harmonics. Comparing the two algorithms using these 
tables, we can conclude that for Case 3 the ALP algorithm gives smaller magnitude es- 
timation errors in most instances than the LMS algorithm for both angles. 
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AVERAGE ESTIMATION ERROR C DEG ) 

-45.0 -3S.0 -27.0 -18.0 -9.0 0.0 9.0 18.0 27.0 3S.0 4S.0 




CASE: LMS3 

EREO: 1000.0 HZ ti- 11 N- 1 1 NITER- 100 

Figure 50. Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 1000 Hz for Case LMS3. 
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AVERnGE ESTIMATION ERROR (OEG) 

-45.0 -35.0 -27.0 -18.0 -9.0 0.0 9.0 18.0 27.0 25.0 4S.0 




COSE: LMS3 

EREQ: 3000.0 HZ H- 1 1 N- 1 1 NITER- 100 



Figure 51. Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 3000 Hz for Case LMS3. 
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Figure 52. Estimates of Direction Cosines UO and VO, Case ALP3: SNR = -3dB and 
q=l 
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Figure 53. Estimates of Direction Cosines UO and VO, Case ALP3: SNR==-3dB and 
q = 3 
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CASE : ALP3 SNR- -3.0 D0 

FREO: 6000.0 HZ H- 1 1 N- 1 1 [-1 P- 10 NITER- 1000 



Figure 54. Estimates of Direction Cosines UO and VO, Case ALP3: SNR = -3dB and 
q = 6 
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Figure 55. Estimates of Direction Cosines UU and VO, Case ALP3: SNR = -6dB and 
q=l 
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Figure 56. Estimates of Direction Cosines U() and VO, Case ALP3: SNR = -6dB and 
q = 3 
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Figure 57. Estimates of Direction Cosines UO and VO, Case ALP3: SNR = -6dB and 
q=6 
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Figure 58. Estimates of Direction Cosines UO and VO, Case ALP3: SNR = -9dB and 
q=l 
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Figure 59. Estimates of Direction Cosines UU and VO, Case ALP3: SNR = -9dB and 
q = 3 
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Figure 60. Estimates of Direction Cosines UO and VO, Case ALP3: SNR = -9dB and 
q = 6 
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AVERAGE ESTIMATION ERROR (DEG) 

-45.0 -36.0 -27.0 -18.0 -9.0 0.0 9.0 18.0 27.0 36.0 45.0 




CASE: ALP3 

FREO: 1000.0 HZ h- 1 1 N- 1 1 I- 1 P- 10 NITER- 1000 

Figure 61. Average Estimation Errors of Depression and Bearing Angles vs. SNR 

at 1000 Hz for Case ALP3. 
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AVERAGE ESTIMATION ERROR (DEG) 

- 45.O - 35.0 - 27.0 - 18.0 - 9.0 0.0 9.0 18.0 27.0 36.0 45.0 




CASE: ALP3 

FREO: 3000.0 HZ M- 1 1 N- 1 1 I- 1 P- 10 NITER- 1000 

Figure 62. Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 3000 Hz for Case ALP3. 



87 



AVERAGE ESTIMATION ERROR (DEG) 

- 45.0 - 36.0 - 27.0 - 18.0 - 9.0 0.0 9.0 18.0 27.0 36.0 45.0 




CASE : ALP3 

FREO: 6000.0 HZ M- 1 1 N- 1 1 I- 1 P-10 NITER- 1000 

Figure 63. Average Estimation Errors of Depression and Bearing Angles vs. SNR 
at 6000 Hz for Case ALP3. 
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D. CASE 4 



Case 4 corresponds to three broadband targets being present two of which have a 
common harmonic. Target 1 was located at (6 0 = 45°, t// 0 = 0°). Target 2 was located at 
(0 O = 45°, i/'o = ISO 0 ). Target 3 was located at ( 0 O =33° , ^ 0 = 47°). A total of six har- 
monics were present in the output electrical signal. The following three subcases were 
examined: 

• Case 4-A: Three targets, two of which are located in the same plane. Targets 1 and 
2 have a common spectral line at 1000 Hz. 

• Case 4-B: Three targets, two of which are located in the same plane. Targets 1 and 
2 have a common spectral line at 3000 Hz. 

• Case 4-C: Three targets two of which are located in the same plane. Targets 1 and 
2 have a common spectral line at 6000 Hz. 

For these cases, only no noise results were compiled for both algorithms in order to 
examine the ability of target discrimination when two targets are in the same plane and 
have a common harmonic. 

Figures 64, 65 and 66 present estimates of the direction cosines U 0 and V 0 for each 
common harmonic, that is at 1000 Hz for Case ALP4-A, 3000 Hz for Case ALP4-B and 
at 6000 Hz for Case ALP4-C. Clearly, we can observe that except for case ALP4-A, the 
algorithm can discriminate the location of two targets. 

Tables 7 through 12 in the Appendix A present numerical data for all the subcases 
of Case 4 for both algorithms. From these tables we can observe that for all the har- 
monics, except the common one, both algorithms correctly identified the location of the 
targets with zero degrees estimation errror. However, for the common harmonic, the 
LMS algorithm located a false fourth target, which was located exactly between the ac- 
tual locations of targets 1 and 2. On the other hand the ALP algorithm, except for the 
subcase ALP4-A, for the common harmonic, identifies two different estimates of the di- 
rection cosine L' 0 from the transfer function plots, which results in two different target 
locations. 
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Figure 64. Estimates of Direction Cosines UO and VO, Case ALP4-A: No Noise and 
q=l 



90 




u 




FREO: 3000.0 HZ H- 11 N- 1 1 1-1 P- 10 NITER- 1000 

Figure 65. Estimates of Direction Cosines UU and VO, Case ALP4-B: No Noise and 
q=3 
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Figure 66. Estimates of Direction Cosines UO and VO, Case ALP4-C: No Noise and 
q = 6 
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E. CASE 5 



Case 5 placed three broadband targets at random locations. Target 1 was located 
at (0 Q = 49°, = 38°). Target 2 was located at (0 O = 5°, i j/ 0 = 137°). Target 3 was located 

at (0„ = 77°, \jj g = 307°). A total of six harmonics were present in the output electrical 
signal. The following three subcases were examined: 

• Case 5-A: Three targets located at random positions. Targets 1 and 2 have a 
common spectral line at 1000 Hz. 

• Case 5-B: Three targets located at random positions. Targets 1 and 2 share a 
common spectral line at 3000 Hz. 

• Case 5-C: Three targets located at random positions. Targets 1 and 2 share a 
common spectral line at 6000 Hz. 

For these cases, only no noise results were compiled for both algorithms in order to 
examine the ability of target discrimination when two targets are in random positions 
and have a common harmonic. 

Figures 67, 68 and 69 present estimates of the direction cosines L’ 0 and V g for each 
common harmonic, that is. at 1000 Hz for Case ALP5-A, at 3000 Hz for Case ALP5-B 
and at 6000 Hz for case ALP5-C. Clearly, we can observe that exept for case ALP5-A, 
the algorithm can discriminate the location of two targets. 

Tables 13 through 18 in the Appendix A present numerical data for all the subcases 
of Case 5 for both algorithms. From these tables we can observe that for all the har- 
monics, except the common one, both algorithms correctly identified the location of the 
targets with zero degrees estimation errror. However, for the common harmonic, the 
LMS algorithm located a false fourth target. On the other hand the ALP algorithm, 
except for subcase ALP5-A, for the common harmonic, identifies two different estimates 
of the direction cosines U 0 and F 0 from the transfer function plots, which results in two 
different target locations. 
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Figure 67. Estimates of Direction Cosines UO and VO, Case ALP5-A: No Noise and 
q= 1 



94 




u 





















A 


\ 


\ 






a 


\ 
































\ 




y 




























































































































































































































































































* 



















































1. 0-0.9 -0.8 -0.7 -0.0 -o'.S -0.1 -0.3 -0.2 -o.l 0.0 0.1 0.2 0.3 0.1 0.5 0.6 0.7 0.8 0.9 1.0 



CR5E:flLP5-0 NO N0I5E 

FREO: 3000.0 HZ fi- 1 1 N- II l- l P’ 10 NITER- 1000 

Figure 68. Estimates of Direction Cosines UO and VO, Case ALP5-B: No Noise and 
q=3 
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Figure 69. Estimates of Direction Cosines U() and VU, Case ALP5-C: No Noise and 
q = 6 
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IV. CONCLUSIONS AND RECOMMENDATIONS 



The goal of this thesis was to evaluate the multiple broadband target localization 
capability and the full angular coverage capability of two frequency domain adaptive 
beamforming algorithms. The main advantage of both algorithms is that no a priori 
knowledge is assumed about the statistics of the received signal, the total number of 
targets present, where they are located, and the frequency content of the targets. 

The two frequency domain adaptive beamforming algorithms evaluated were the 
modified complex Least-Mean-Square (LMS) algorithm developed by Ziomek and Chan 
[Ref. 6], and the modified complex Adaptive Linear Prediction-Error (ALP) algorithm 
which was presented for first time. Chapter 3 presented the results of five test cases 
which were designed to test both algorithms for several different capabilities. It should 
also be noted here that the number of time samples taken per element of the array never 
exceeded 65 time samples, an amount which represents a small number of data points. 

Summarizing the simulation results discussed in Chapter 3, we can list the following 
conclusions: 

• In most cases tested, both algorithms can localize multiple broadband targets and 
they have full angular coverage capability. 

• Increasing either the sampling rate or the array size yields better estimates in lower 
SNR for both algorithms 

• ALP generally gives better estimates in lower SNR than LMS 

• ALP can, in most of the "no noise" cases, discriminate between two targets that 
have a common harmonic while LMS cannot 

• ALP is more computationally intensive than LMS and requires much more iter- 
ations to approach a steady-state error. 

• ALP is more flexible than LMS because it can process only a part of the total 
number of array elements to localize the target 

In the course of this investigation, several possible areas for future research pre- 
sented themselves: 
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• further study of the effects of varying the step-size parameter, n , in an attempt to 
decrease the magnitude of the estimation error. 

• investigation of parallel processing implementation for the ALP algorithm in order 
to reduce the computation time. 

• investigation of larger array size in lower SNR mainly using the ALP algorithm 

• investigation of a noise reduction system prior to processing the output electrical 
signals. 
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APPENDIX A. 



This appendix presents the numerical data for all of the test cases. Tables 1 through 
6 present numerical data for the three first "noisy" cases, Cases 1 through 3, for both 
algorithms, while tables 7 through 12 present numerical data from the two "no noise" 
cases, Case 4 and 5, for both algorithms. The following notation is used in all the tables: 

• LMS: When the frequency domain modified complex Least-Mean-Square adaptive 
beamforming algorithm is used. 

• ALP: When the frequency domain modified complex adaptive linear prediction- 
error filter beamforming algorithm is used 

• 6 0 : Depression Angle 

• >j > 0 : Bearing Angle 

• c 9q : Average Estimation Error for Depression Angle 

• e v>0 : Average Estimation Error for Bearing Angle 
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Table 1. NUMERICAL DATA CORRESPONDING TO CASE LMS1 



Case 1 : One target located at broadside relative to planar array - six harmonics 
present in output electrical signal - LMS algorithm is used - 100 iterations 





Location (deg) 


Average Estimation Error (deg) 


Harmonic (q) 






SNR 


0o 


<A o 


-3 dB 


-6 dB 


-9 dB 














1 


0 


270 


-2.169 


-1.7S9 


-4.382 


2 


0 


270 


-0.S10 


-1.426 


-2.170 


3 


0 


270 


-0.561 


-0.607 


-1.233 


4 


0 


270 


-0.360 


-0.760 


-0.649 


5 


0 


270 


-0.362 


-0.363 


-0.848 


6 


0 


270 


-0.299 


-0.435 


-1.596 



Table 2. NUMERICAL DATA CORRESPONDING TO CASE ALP1 



Case 1: One target located at broadside relative to planar array - six harmonics 
present in output electrical signal - ALP algorithm is used - 1000 iterations 





Location (deg) 


Average Estimation Error (deg) 


Harmonic (q) 






SNR 


00 


<Ao 


-3 dB 


-6 dB 


-9 dB 












g »p 


1 


0 


270 


-3.623 


-4.444 


-6.056 


2 


0 


270 


-1.307 


-2.739 


-2.404 


3 


0 


270 


-1.325 


-1.831 


-1.993 


4 


0 


270 


-1.231 


-1.429 


-1.463 


5 


0 


270 


-0.753 


-0.740 


-1.746 


6 


0 


270 


-0.471 


-0.665 


-0.951 



100 



Table 3. NUMERICAL DATA CORRESPONDING TO CASE LMS2 



Case 2: Three targets located at random positions - six harmonics present in output 
electrical signal - LMS algorithm is used - 100 iterations 



Harmonic 

(q) 


Location (deg) 


Average Estimation Error (deg) 


0o 


<Ao 


SNR 


-3 dB 


-6 dB 


-9 dB 


% 


e 

v 'n 








A'n 


i 


49 


38 


0.265 


-0.124 


-0.411 


-0.832 


-0.571 


0.325 


2 


5 


137 


0.077 


-2.714 


-0.161 


-5.017 


0.097 


3.451 




77 


307 


3.037 


-0.108 


0.351 


-2.625 


7.991 


9.333 


4 


77 


307 


-10.532 


1.677 


-5.032 


0.019 


3.858 


26.436 


5 


5 


137 


-0.022 


0.560 


-0.155 


1.376 


0.07S 


0.002 


6 


49 


38 


1.0S1 


-1.531 


-10.852 


4.382 


-2.355 


-8.7S4 



Table 4. NUMERICAL DATA CORRESPONDING TO CASE ALP2 



Case 2: Three targets located at random positions - six harmonics present in output 
electrical signal - ALP algorithm is used - 1000 iterations 



Harmonic 

(q) 


Location (deg) 


Average Estimation Error (deg) 


e 0 


1 0 


SNR 


-3 dB 


-6 dB 


-9 dB 


^0 


... ^ ... 




e *n 






1 


49 


38 


-1.470 


0.098 


-0.826 


0.3S0 


8.243 


-2.443 


2 


5 


137 


0.134 


6.094 


-1.091 


-8.072 


0.S79 


1.252 


3 


77 


307 


-1.72S 


0.082 


-3.247 


-0.361 


-0.697 


1.131 


4 


77 


307 


-1.968 


0.526 


-0.964 


0.218 


-1.763 


-0.035 


5 


5 


137 


0.019 


-0.191 


0.180 


1.267 


0.830 


-6.591 


6 


49 


38 


-0.173 


0.146 


-0.197 


-0.019 


2.06S 


-5.109 
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Table 5. NUMERICAL DATA CORRESPONDING TO CASE LMS3 



Case 2: One target located at endfire relative to planar array - six harmonics present 
in output electrical signal - LMS algorithm is used - 100 iterations 



Harmonic 

(q) 


Location (deg) 


Average Estimation Error (deg) 


00 


'I'o 


SNR 


-3 dB 


-6 dB 


-9 dB 


e e n 




e #n 


% 






1 


90 


90 


3.355 


-0.055 


20.000 


0.545 


21.186 


1.023 


2 


90 


90 


33.033 


-0.200 


40.425 


-0.646 


27.226 


-0.330 


J 


90 


90 


19.490 


-0.019 


16.401 


-0.229 


16.706 


-0.325 


4 


90 


90 


2.247 


0.006 


7.365 


0.287 


12.025 


-0.139 


5 


90 


90 


2.136 


-0.0S3 


1.146 


0.098 


14.692 


-18.528 


6 


90 


90 


70.S77 


-84.120 


69.426 


-64.811 


46.737 


-90.393 



Table 6. NUMERICAL DATA CORRESPONDING TO CASE ALP3 



Case 2: One target located at endfire relative to planar array - six harmonics present 
in output electrical signal - ALP algorithm is used - 1000 iterations 



Harmonic 

(q) 


Location (deg) 


Average Estimation Error (deg) 


00 


»Ao 


SNR 


-3 dB 


-6 dB 


-9 dB 




e *n 




e *o 


€ff n 


e *‘fl 


i 


90 


90 


3.948 


0.217 


7.806 


-0.136 


6.672 


-0.984 


2 


90 


90 


7.287 


-0.362 


11.323 


-0.4S8 


9.034 


0.329 


3 


90 


90 


5.342 


-0.107 


3.732 


-0.370 


1.935 


-0.767 


4 


90 


90 


2.233 


-0.113 


2.145 


0.286 


7.725 


-0.023 


5 


90 


90 


0.986 


-0.160 


3.512 


0.076 


14.090 


-6.512 


6 


90 


90 


5.229 


-89.997 


7.051 


-90.029 


16.548 


-89.0S8 
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Table 7. NUMERICAL DATA CORRESPONDING TO CASE LMS4-A 



Case 4a: Three targets two of which are located in the same plane - targets one and 
two have a common spectral line at 1000 Hz - six harmonics present in output elec- 
trical signal - LMS algorithm is used - 100 iterations - no noise 


Harmonic (q) 


Location (deg) 


Average Estimation Error 
(deg) 


00 




e Jn 




1 


45 


0 


45.000 


-90.000 


1 


45 


180 


45.000 


90.000 


2 


45 


0 


0.000 


0.000 


3 


33 


47 


0.000 


0.000 


4 


33 


47 


0.000 


0.000 


5 


45 


180 


0.000 


0.000 


6 


33 


47 


0.000 


0.000 



Table 8. NUMERICAL DATA CORRESPONDING TO CASE ALP4-A 



Case 4a: Three targets two of which are located in the same plane - targets one and 
two have a common spectral line at 1000 Hz - six harmonics present in output elec- 
trical signal - ALP algorithm is used - 1000 iterations - no noise 


Harmonic (q) 


Location (deg) 


Average Estimation Error 
(deg) 


00 








1 


45 


0 


-44.921 


0.0 


1 


45 


180 


-44.921 


180.000 


2 


45 


0 


0.009 


0.000 


3 


33 


47 


0.037 


-0.011 


4 


33 


47 


0.037 


-0.011 


5 


45 


180 


0.009 


0.000 


6 


33 


47 


0.037 


-0.011 
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Table 9. NUMERICAL DATA CORRESPONDING TO CASE LMS4-B 



Case 4b: Three targets two of which are located in the same plane - targets one and 
two have a common spectral line at 3000 Hz - six harmonics present in output elec- 
trical signal - LMS algorithm is used - 100 iterations - no noise 


Harmonic (q) 


Location (deg) 


Average Estimation Error 
(deg) 


#0 


< A 0 




e *'f\ 


1 


33 


47 


0.000 


0.000 


2 


45 


0 


0.000 


0.000 


3 


45 


0 


45.00 


-90.00 


3 


45 


180 


45.0 


90.0 


4 


33 


47 


0.000 


0.000 


5 


45 


ISO 


0.000 


0.000 


6 


33 


47 


0.000 


0.000 



Table 10. NUMERICAL DATA CORRESPONDING TO CASE ALP4-B 



Case 4b: Three targets two of which are located in the same plane - targets one and 
two have a common spectral line at 3000 Hz - six harmonics present in output elec- 
trical signal - ALP algorithm is used - 1000 iterations - no noise 


Harmonic (q) 


Location (deg) 


Estimation Error (deg) 


00 


^0 


.... %. 




1 


33 


47 


0.037 


-0.011 


2 


45 


0 


0.009 


0.000 


3 


45 


0 


1.56 


0.000 




45 


180 


1.56 


0.000 


4 


33 


47 


0.032 


-0.011 


5 


45 


180 


0.009 


0.000 


6 


33 


47 


0.037 


-0.011 
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Table 11. NUMERICAL DATA CORRESPONDING TO CASE LMS4-C 



Case 4c: Three targets two of which are located in the same plane - targets one and 
two have a common spectral line at 6000 Hz - six harmonics present in output elec- 
trical signal - LMS algorithm is used - 100 iterations - no noise 


Harmonic (q) 


Location (deg) 


Average Estimation Error 
(deg) 


00 


•Ao 




^'n 


1 


33 


47 


0.000 


0.000 


2 


45 


0 


0.000 


0.000 


3 


33 


47 


0.000 


0.000 


4 


33 


47 


0.000 


0.000 


5 


45 


ISO 


0.000 


0.000 


6 


45 


0 


45.000 


-90.000 


6 


45 


ISO 


45.000 


90.000 



Table 12. NUMERICAL DATA CORRESPONDING TO CASE ALP4-C 



Case 4c: Three targets two of which are located in the same plane - targets one and 
two have a common spectral line at 6000 Hz - six harmonics present in output elec- 
trical signal - ALP algorithm is used - 1000 iterations - no noise 


Harmonic (q) 


Location (deg) 


Average Estimation Error 
(deg) 




<Ao 


% 


_ - - 


1 


33 


47 


0.037 


-0.011 


2 


45 


0 


0.009 


0.000 


3 


33 


47 


0.037 


-0.011 


4 


33 


47 


0.037 


-0.011 


5 


45 


180 


0.009 


0.000 


6 


45 


0 


2.070 


0.000 


6 


45 


ISO 


2.070 


0.000 
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Table 13. NUMERICAL DATA CORRESPONDING TO CASE LMS5-A 



Case 5a: Three targets located at random positions - targets one and two have a 
common spectral line at 1000 Hz - six harmonics present in output electrical signal - 
LMS algorithm is used - 100 iterations - no noise 


Harmonic (q) 


Location (deg) 


Average Estimation Error 
(deg) 


e o 


'/'o 


% 


e *n 


1 


49 


38 


-13.000 


252.594 


1 


5 


137 


-41.000 


-16.406 


2 


5 


137 


0.000 


0.000 


J 


77 


307 


0.000 


0.000 


4 


77 


307 


0.000 


0.000 


5 


5 


137 


0.000 


0.000 


6 


49 


38 


0.000 


0.000 



Table 14. NUMERICAL DATA CORRESPONDING TO CASE ALP5-A 



Case 5a: Three targets located at random positions - targets one and two have a 
common spectral line at 1000 Hz - six harmonics present in output electrical signal - 
ALP algorithm is used - 1000 iterations - no noise 


Harmonic (q) 


Location (deg) 


Average Estimation Error 
(deg) 


00 


lAo 






1 


49 


38 


27.121 


-6.674 


1 


5 


137 


-16.879 


92.326 


2 


5 


137 


-0.038 


-o.oos 


3 


77 


307 


0.094 


0.012 


4 


77 


307 


0.094 


0.012 


5 


5 


137 


0.006 


-0.328 


6 


49 


38 


0.094 


0.012 
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Table 15. NUMERICAL DATA CORRESPONDING TO CASE LMS5-B 



Case 5b: Three targets located at random positions - targets one and two have a 
common spectral line at 3000 Hz - six harmonics present in output electrical signal - 
LMS algorithm is used - 100 iterations - no noise 


Harmonic (q) 


Location (deg) 


Average Estimation Error 
(deg) 


Op 


4> o 


e Q 


% 


1 


77 


307 


0.000 


0.000 


2 


5 


137 


0.000 


0.000 


3 


49 


38 


-13.000 


252.594 


3 


5 


137 


-41.000 


-16.406 


4 


77 


307 


0.000 


0.000 


5 


5 


137 


0.000 


0.000 


6 


49 


38 


0.000 


0.000 



Table 16. NUMERICAL DATA CORRESPONDING TO CASE ALP5-B 



Case 5b: Three targets located at random positions - 

targets one and two have a common spectral line at 3000 Hz - six harmonics present 
in output electrical signal - ALP algorithm is used - 1000 iterations - no noise 


Harmonic (q) 


Location (deg) 


Estimation Error (deg) 




’/'o 






1 


77 


307 


0.094 


0.012 


2 


5 


137 


-0.038 


-0.008 


3 


49 


3S 


2.94 


0.39 


3 


5 


137 


1.505 


-1.992 


4 


77 


307 


0.094 


0.012 


5 


5 


137 


0.006 


-0.328 


6 


49 


38 


0.094 


0.012 
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Table 17. NUMERICAL DATA CORRESPONDING TO CASE LMS5-C 



Case 5c: Three targets located at random positions - targets one and two have a 
common spectral line at 6000 Hz - six harmonics present in output electrical signal - 
LMS algorithm is used - 100 iterations - no noise 


Harmonic (q) 


Location (deg) 


Average Estimation Error 
(deg) 


0, 


'Ao 




% 


1 


77 


307 


0.000 


0.000 


2 


5 


137 


0.000 


0.000 


3 


49 


3S 


0.000 


0.000 


4 


77 


307 


0.000 


0.000 


5 


5 


137 


0.000 


0.000 


6 


5 


137 


-41.000 


-16.406 


6 


49 


38 


-13.000 


252.594 



Table 18. NUMERICAL DATA CORRESPONDING TO CASE ALP5-C 



Case 5c: Three targets located at random positions - targets one and two have a 
common spectral line at 6000 Hz - six harmonics present in output electrical signal - 
ALP algorithm is used - 1000 iterations - no noise 


Harmonic (q) 


Location (deg) 


Average Estimation Error 
(deg) 


0o 


<Ao 


e fi 




1 


77 


307 


0.094 


0.012 


2 


5 


137 


-0.03S 


-0.008 


3 


49 


3S 


0.094 


0.012 


4 


77 


307 


0.094 


0.012 


5 


5 


137 


0.006 


-0.328 


6 


5 


137 


2.615 


-3.84 


6 


49 


3S 


0.49S 


2.516 
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